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We study the possibility that large scale cosmological perturbations have been generated during the domination 
and decay of a massive particle species whose mass depends on the expectation value of a light scalar field. We 
discuss the constraints that must be imposed on the field in order to remain light and on the annihilation cross 
section and decay rate of the massive particles in order for the mechanism to be efficient. We compute the resulting 
curvature perturbations after the mass domination, recovering the results of Dvali, Gruzinov, and Zaldarriaga in 
the limit of total domination. By comparing the amplitude of perturbations generated by the mass domination to 
those originally present from inflation, we conclude that this mechanism can be the primary source of perturbations 
only if inflation does not rely on slow-roll conditions. 


1. Introduction 

Although inflation has recently become the 
paradigm of the early universe and the domi¬ 
nant contender for generating the observed large 
scale density perturbations pQ, it is known to be 
plagued by the difficulty of being embedded in 
most current particle physics models. 

The usual hypothesis is that primordial density 
perturbations originate from the vacuum quan¬ 
tum fluctuations of the inflaton field. The fact 
that the inflaton has to remain light during infla¬ 
tion and that its vacuum energy is constrained 
by the amplitude of the observed cosmological 
perturbations make it difficult to construct moti¬ 
vated models |21 ■ For this reason, scenarios where 
inflation does not rely on slow-roll conditions have 
recently been proposed (see e.g., PI). 

If, as in these scenarios, the inflaton has no fluc¬ 
tuations, it is necessary to find a mechanism to 
generates cosmological perturbations. In the cur- 
vaton scenario m curvature perturbations are 
generated by the late decay of a light scalar field 
whose fluctuation S(f> is amplified during in¬ 
flation with a quasi-scale invariant spectrum, and 
amplitude set by the Hubble rate at horizon cross¬ 
ing, Scj) ^ H. This mechanism opens up new pos¬ 
sibilities for observations and for model building. 
In particular, it can liberate the inflaton from the 
requirement that it is responsible for the primor¬ 
dial perturbations. 


In the same spirit, inspired by the idea that 
coupling constants and masses of particles dur¬ 
ing the early universe may depend on the value 
of some light scalar field, a number of people 
have proposed that perturbations could be gener¬ 
ated from the fluctuations of the inflaton coupling 
to ordinary matter during reheating, a mecha¬ 
nism called inhomogeneous reheating or modu¬ 
lated fluctuations [HI (see also PHI and piT)] for 
the perturbations). 

Here we study the possibility that cosmological 
perturbations are generated during the phase of 
domination of a massive particle species '0 whose 
mass is modulated by the value of a light field 
0, as originally proposed in m These parti¬ 
cles are produced and thermalize after inflation. 
If sufficiently long-lived, they can freeze out and 
eventually dominate the universe when they be¬ 
come non-relativistic. Due to the mass fluctua¬ 
tion, the mass domination process becomes inho¬ 
mogeneous and converts the fluctuation of 0 into 
curvature perturbations. 

In we studied the perturbations generated 
by this mechanism making use of the formalism 
introduced in m for interacting fluids. Here we 
take a different approach and, in particular, we 
consider more appropriate initial conditions for 
the perturbations of the massive particles after 
inflation. Due to this different choice of initial 
conditions, the result found here for the final cur¬ 
vature perturbation in terms of the mass fluctua- 
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tions does not agree with m, but confirms m 
In the next section we model the inhomoge¬ 
neous mass domination. We derive and discuss 
the conditions that must be imposed on the field 
for remaining light (Sec. |2l and on the annihila¬ 
tion cross section and decay rate of the massive 
particles in order for tjj to dominate the universe 
and for the mechanism to be efficient (Sec. @J. 
Finally, we study the perturbations and the ob¬ 
servational consequences of this mechanism, and 
we discuss when it can successfully liberate infla¬ 
tion. 

2. Modeling the mass fluctuation 

In this section we discuss the homogeneous 
evolution equations governing a fluid of non- 
relativistic particles ip with a rest mass that de¬ 
pends on the value of a scalar field. Since we want 
to describe the universe after reheating, we also 
add a radiation fluid and let ■i/' to decay into radi¬ 
ation. The derivation given here of the coupling 
between ip and (/) is ‘intuitive’: A more rigorous 
treatment can be found in the literature on scalar- 
tensor theories such as in mm- 
We consider a flat Friedmann-Lemaitre- 
Robertson-Walker universe, with metric ds^ = 
—dt^ + a?{t)dx^. The Hubble rate H = a/a 
{' = d/dt) is given by 

= 2^2 + ^7 + P<t>)^ ( 1 ) 

where px is the energy density of the X species 
and mp = (SttG)”^/^ is the reduced Planck mass. 

The non-relativistic particles considered here 
are similar to cold dark matter particles: they 
are collisionless and non-relativistic. Hence their 
pressure vanishes and their energy density is given 
by the product of the particle number density and 
the mass, 

, ( 2 ) 

where u^, is the four-velocity. The mass of tp de¬ 
pends on a scalar field </, 

We consider ip after freeze-out so that np is 
not changed by annihilation processes. However, 
since we want to allow -ip to decay into radiation, 
we define a decay rate of change of the number 


of Ip particles as V^(n,^u^) = —Tn^. With this 
equation, Eq. m yields a conservation equation 
for the energy of ip, written in covariant form na, 

= (3) 

where we have defined the (dimensionless) mass 
coupling function a = d\nmp/dlii(p. Even in the 
absence of coupling to radiation, due to the inter¬ 
action with cp, the energy of ip is not conserved. 

Since ip decays only into radiation, the evolu¬ 
tion equation governing the radiation energy den¬ 
sity is 

= -Tpp. (4) 

The background conservation equations for the 
two fluids then become (with P-y = p^/i) 

Pp = -{iH - a^/(p + T)pp, (5) 

p^ = -AHp^ + Tpp. (6) 

The scalar field (p has a standard kinetic term 
and a potential V{(p). However, its evolution 
equation is sourced by an additional term due to 
the coupling to the bath of non-relativistic parti¬ 
cles Ip. Indeed, the requirement that the sum of 
the three energy-momentum tensors is conserved. 
Ex = 0, yields 

X^d>^iP = V' + {a/cP)pp. ( 7 ) 

The background Klein-Gordon equation for the 
scalar field in discussed in the following section. 

3. Scalar field behavior 

The background part of Eq. o can be rewrit¬ 
ten as 

P + 3HP + V' =-{a/(p)pp, ( 8 ) 

where ' = d/d(p. Taking into account the coupling 
to Ip, the scalar field behaves as a field in the effec¬ 
tive potential Veff{(p) = V{(p)+pp. The coupling 
can have the effect of increasing the mass of (p 
and drive its evolution. 

Exploring the behavior of the scalar field dur¬ 
ing the phase that goes from the radiation domi¬ 
nated era to the mass domination requires a form 
for V{(p) and mp{(p). Several form of these func¬ 
tions have been explored in the literature. If the 
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effective mass of (j) is made larger than H by the 
coupling to V') the field rapidly sets to the mini¬ 
mum of VeS mm- Here, however, we want to 
protect the mass of (j) and maintain the field light, 
to conserve the primordial fluctuations inherited 
from inflation CH. This puts strong constraints 
on the initial value of (f). 

We consider a massive field, V{(j)) = , 

with 7710 < r, so that < H until the decay of 
Ip. The condition for the field to remain light is 
then 

(9) 

where is the abundance of ip, = pv,/p, at 
its decay. A more tighter condition is to impose 
that V' always dominates over {a/(p)p^, 

(pin/mp > ( 10 ) 

If this condition holds at some time ‘in’ deep in 
the radiation era, it is maintained throughout the 
domination. These conditions can be realized for 
(p ~ TTip if is small enough. An analogous 
condition on (p was first advocated in CH in or¬ 
der to protect its mass from acquiring too large 
thermal corrections due to the coupling with the 
bath of Ip particles. 

If, on the other hand, (p is too large, the scalar 
field can dominate the universe by entering a sec¬ 
ond inflationary era. We must hence require that 

(p/mp < r/?7i0. (11) 

In the following we assume the stronger condi¬ 
tion cni). This simplifies considerably the treat¬ 
ment of perturbations than by only requiring the 
weaker condition, Eq. ®. One can indeed study 
the background evolution of (p neglecting the cou¬ 
pling to Ip. Equation © becomes 

^ 0, a(t)(xt^. ( 12 ) 

The solution of this equation (regular for m^t —>■ 
0) is given in terms of the Bessel function Ji, with 
77 = (3/3— 1)/2: (p ~ (p^A„J^{m^t)/{m^tY, where 
Ai, is a numerical constant and the star stands for 
the value during inflation. Hence we obtain 

3/3 

3H^- < 1. (13) 


Note that this solution is not of ‘slow-roll type’: 
The acceleration (p is not small with respect to 
H(p. However, as during slow-roll, the kinetic en¬ 
ergy of the scalar field is subdominant with re¬ 
spect to the potential energy, ~ —Pij), while 
the adiabatic speed of sound depends on /3, ~ 

-l-2/(3/3). 

4. Conditions for the mass domination 

For the mass fluctuations to be imprinted into 
the density perturbations it is essential that ip 
dominates or at least becomes a significant com¬ 
ponent of the energy density of the universe be¬ 
fore its decay. If ip froze out at a temperature T 
such that m.p/T was not much larger than 1, then 
the species can have a significant relic abundance 
and, if sufficiently long-lived, can eventually dom¬ 
inate the universe before decaying. Here we want 
to study when this condition is realized. 

As shown by Eq. m, the evolution of the light 
field is slow with respect to the expansion rate, 
—cp/(p {m^f,/HY H ^ H, so that the mass vari¬ 
ation is also slow, 7770/7710 ^ H. Furthermore, 
we are considering non-relativistic particles, i.e., 
T < m.,p ^ H < 777 ^/ 777 p ^ 7770 . Thus, when 
studying their abundance, the particles ip can be 
considered as having constant mass m, and we 
can consistently neglect the field coupling on the 
right hand side of Eq. ©• 

The system of equations © and © then de¬ 
scribes the decay of a dust fluid into radiation and 
can be easily solved numerically (see e.g., mm)- 
The initial conditions for this system are taken 
well into the radiation era, H 0 ^ I, once ip has 
frozen out. As shown in Ref. close to the 
initial condition H 0 <C 1, 120 — p(r/iJ)^/^. The 
value of the parameter 

p = (14) 

determines the trajectories followed by the system 
on the ( 00 , H)-plane. In particular it determines 
the abundance of ip at its decay, O^®®. This is 
shown in Fig. ^ as a function of p. If initially 
r 3?0|, i.e., p <C I, then the decay is almost 
instantaneous and ip does not have the time to 
dominate, O^®® <C 1. On the contrary, when p is 
initially large, then O^®® = 1. Thus, the condition 
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Figure 1. The dashed line shows the abundance 
of Ip at its decay, as a function of the initial 
parameter p. This is compared to the efficiency 
parameter r (solid line), defined in Sec. 15.21 


for Ip to dominate the universe becomes p ^ 1 . 

Now we want to translate this into a condition 
on the decay rate and annihilation cross section 
of pj. On using the definition of p, Eq. m , we 
can write the condition for the domination as 


P- g 


1 /4 


•s/mpT 


> 1 , 


(15) 


where g is the number of relativistic degrees of 
freedom and n^/s is the relic abundance of the tp 
particles when they freeze-out - s is the entropy 
density. In order to derive m we have used p ~ 
p^ ~ gT'^, and s ~ gT^- 

After freeze-out and for T H, the relic abun¬ 
dance n^/s is constant (we assume g constant). 
The relic abundance can be evaluated by solving 
the Boltzmann equation describing the freeze out 
of the annihilation of ip particles and antiparticles 
with cross section a a, which we take to be inde¬ 
pendent of the particles energy. In this case there 
is an approximate solution for the abundance at 
freeze-out im, 


^_ 100 

s mpmpg^/^{crAv)’ 


(16) 


where {aAv) is the thermal average of the total 
cross section times the relative velocity v. On 
using this solution, Eq. US) becomes independent 
of the mass m, 


100 /mpNl/2 TOp^ ^ 

^ g^/'^ V r / {(Tav) 


(17) 


This relation holds if the ^-particles are subdom¬ 
inant at freeze-out. The smaller the annihilation 
cross section is, the earlier ip freeze out before 
becoming non-relativistic and the larger is 

On requiring that the massive particles decay 
before nucleosynthesis, i.e., F > (lMeV)^/mp, we 
find that only for (aAv) < 10 /p^/^( 10 ^GeV)“^ 
does Ip dominate the universe. Hence, the initial 
thermal equilibrium by annihilation of particles 
and antiparticles must be maintained by some 
gauge interaction which is weaker than those of 
the standard model. This excludes that ip is made 
of standard model particles. 


5. Cosmological perturbations from the 
mass domination 

In this section we derive the coupled perturbed 
equations of motion for the two fluids and the 
scalar field. Then we discuss their solutions, an¬ 
alytically and numerically. The aim is to study 
the mechanism of conversion of scalar field fluc¬ 
tuations into curvature perturbations during the 
mass domination. 

We describe scalar perturbations in the metric 
with line element 

ds^ = —(1 -I- 2^)dt^ + 0^(1 — 2^)Sijdx^dx^ , (18) 

which follows from the absence of anisotropic 
stress perturbations. The quantity 4) corresponds 
to the Bardeen potential in longitudinal gauge. 

In order to perturb the energy-momentum ten¬ 
sor of the three components ip, 7 , and (p, we in¬ 
troduce the energy density and pressure pertur¬ 
bations, 5px and 6Px, and the scalar field pertur¬ 
bation 5(p. It is also useful to define the relative 
perturbation 6x = dpx/px- In terms of these 
quantities, we perturb Eqs. 0 , 0 , and o, and 
since we are interested in the large scale pertur¬ 
bations, we drop the gradient terms, that vanish 
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in the large scale limit. This yields 

,5^ - 3$ = - r$, (19) 

j^-4$ = r^(5^-,5^ + $), (20) 

Pi 

54> + "iHdfj) + _ 

4$<^-2y'$-(a/</.)p^(5^+2$). (21) 

The perturbed energy constraint equation reads 

6i/($ + H^) = -mp"^ (pjjjS.^ + pjSj + p^d^).{22) 

Equation ED can be opportunely simplified: 
Since the scalar field is light and subdominant 
we can neglect the right hand side. Assuming 
the strong constraint (Cnj we can also neglect 
the mass correction due to the coupling. Equa¬ 
tion ED then simplifies to the equation of a light 
test field, 

(5> -h SHScj) + ml5(j) = 0, (23) 

that in the massive case is the same as the equa¬ 
tion for the background field. Hence, the ratio 
between 5(j) and (j) is constant in time, 

5(j)/4> = 54)^, / (24) 

Assuming only the weaker constraint implies 
that 5(j)/4> may slowly vary, which makes the 
mechanism more difficult to study, although it 
maintains its qualitative features. 

5.1. Analytic calculation 

It is possible to analytically solve the system 
of equations ira . CT . and in the case 

that the massive particles dominate completely 
the universe before decaying. We can hence con¬ 
sistently assume that we are far from the decay, 
H ^ r, so that the source terms proportional to 
r in Eqs. (imil and can be neglected for the 
moment. Equations and (HOI) then simplify, 

— 3<i> = a5(j)*/4*j (25) 

= 0, (26) 

where we have used Eq. (1^ . 

We first discuss the initial conditions of these 
equations. They are defined in the radiation dom¬ 
inated era, just after the freeze-out of ip. In this 


limit the Bardeen potential $ is constant, and the 
perturbation of the radiation fluid can be found 
from the constraint 123, on using p ~ p^. This 
yields 5^'^ = —2$*, where is the value of $ 
at the beginning of the radiation dominated era 
(e.g., after inflation). This initial perturbation is 
usually dropped and considered negligible in the 
treatment of models where perturbations are pro¬ 
duced after inflation by a light field. Here how¬ 
ever, in the spirit of m, we retain this term since 
it can turn out to be larger than the curvature 
perturbation generated by the mass domination 
(see below). 

The value of the initial condition is crucial 
to correctly determine the final curvature pertur¬ 
bation produced. Unfortunately it is not unam¬ 
biguously defined. Here we set by assuming 
that radiation and ip particles were created by the 
same mechanism (e.g., reheating) with the same 
number density perturbation,^ 



= 0 


z!p!L - "aW 
^b) 4 • 


(27) 


Thus, considering also the mass fluctuation, the 
relative initial density perturbation of tp reads 


M) _ 


^+3 b) 
.b) A + 4 ^^^ ■ 


(28) 


The initial conditions are thus given by 

= aJ4*/(p* - (3/2)$*, = -2$*. (29) 

Equations (ESI and (ESI) can be easily solved 
with these initial conditions, yielding 


<5(5* 9 

= a^-k3$--$*. 


6^ = 4$-6$*, 


(30) 

(31) 


where the terms proportional to are integra¬ 
tion constants. Only the first of these equations 

^In Ref. 1121 the initial condition 5^/3 = i.e., ab¬ 

sence of relative entropy perturbation between "0 and 7 , 
was used. Only if is constant, is this condition the 
same as Eq. EJ. The choice of Eq. 1771 is at the origin 
of the discrepancy between the result of this work and the 
result of Eu on the final total curvature perturbation. 
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is used in the following. At the final stage of the 
mass domination, when the non-relativistic par¬ 
ticles dominate the universe, the constraint 123 
yields = —2$. Combined with Eq. (HU it im¬ 
plies that the final curvature perturbation is 


9 1 ( 50 * 

■" = To*--5“X 


(32) 


This result holds during the matter [ip) domi¬ 
nated era, before the decay of the massive par¬ 
ticles. It is related to the curvature perturbation 
in the subsequent radiation era after ip decay by 
the well-known relation = (10/9)$^ PO- 
Equation (|23 can be expressed in terms of the 
curvature perturbation defined on the uniform 
density hypersurface, ( = —$-I-(5p/[3(/9-|-P)] 
m, which is known to be conserved for adia¬ 
batic perturbations. In the radiation dominated 
era ( = —(3/2)$^, which yields 


< = + 

o m-ifj 


[nf^ = I). 


(33) 


dec 


The coefficient 1/3 in this equation confirms the 
result of [TT|. 


5.2. Numerical results 

Here we solve Eqs. llTIIll . (1^ . and nu¬ 
merically, so that we can study the case where 
the massive particles do not dominate the uni¬ 
verse completely before decaying. We initially 
consider the particular case of a = constant (e.g., 
= A0"). In this case the first term on the 
right hand side of Eq. null vanishes. The prob¬ 
lem then reduces to the study of a pressureless 
fluid with non-vanishing initial perturbation 123 
that decays into a relativistic fluid 7. 

This problem has been studied in m in the 
case of the curvaton field. There, it was found 
that the final curvature perturbation is C = 

(<I>* = 0 was assumed), where r is an efficiency 
parameter which can be computed as a function 
of p. We have solved Eqs. m and numeri¬ 
cally and the solution is 

C = C* + ^Q!*^, {m^ = X(p, (34) 

with r < 1 represented in Fig. ^ Here an ini¬ 
tial non-vanishing perturbation 0* has also been 


included. Fig. ^ shows that r ~ As ex¬ 

pected, the more ip dominates the universe before 
decaying, the more the mechanism of generation 
of curvature perturbations is efficient. 

If a varies slowly during the phase of domina¬ 
tion and decaying, we can extend the results of 
the numerical calculation to a general a. We find 


r (50* . r 

C — C* + gadec-— — C* + g - 

O O TTlip 


dec 


(35) 


This is our main result. It generalizes the results 
of in to the case where the non-relativistic par¬ 
ticles do not dominate the universe completely 
before decaying into radiation. 


6. Non-Gaussianities 


Here we discuss an unambiguous observational 
consequence of the mass domination, namely the 
presence of non-Gaussianities in the perturba¬ 
tions. The non-Gaussianities generated by the 
mass domination have been first computed in m 
in the case of total domination = 1) (see 

also I2ii2a). The possible presence of isocurva¬ 
ture perturbations has been studied in m 
For convenience we define r = radec — 
H(J,®®Q(dec- We also define the ratio between the 
perturbation generated by the mass domination 
and the primordial curvature perturbation 0* [see 
Eq. 13511 ] . 


1 f 5(p^: 
(p 


We can parameterize the level of non- 
Gaussianities with the non-linear parameter /nl 

defined as C = Cl - (3/5)/nlCl I ^3I24| . where Cl 

represents the Gaussian (linear) contribution to 
the total curvature perturbation C- Motivated by 
inflation, we assume 0* to be Gaussian. 

Non-Gaussianities arise when the field pertur¬ 
bation (50* becomes larger than ^ 1O“^0*. This 
happens when the efficiency is low, f <§C 1. In 
this case, terms quadratic in 50*/0*, that have 
been neglected in the linear calculation, become 
important 


f / (50* 
3 V 0* 



+ 0 *- 


(37) 
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This yields the non-linear parameter 

/NL = -^(l + i?-')-'. (38) 

As expected, when R ^ 0, /nl ^ 0. If, on the 
contrary, R is large, i.e., C is mainly due to the 
mass domination, /nl = ~5/(27) m As no¬ 
ticed in m this number is in the ballpark of 
future experiments like Planck, that will be able 
to detect /nl ^ 5 at 2cr-level j^S]- Therefore 
we may have a well detectable signature of non- 
Gaussianities if r < 0.5. 

Here we are also interested in the lower 
bound on f from the current bound on non- 
Gaussianities. The WMAP limit on /nl corre¬ 
sponds to —58 < /nl < 134 (2(T-level) PK]. which 
translates into f > 0 . 02 . 


7. Liberating inflation? 


As mentioned in the introduction, the main 
reason to use a light field other than the infla- 
ton as a source of curvature perturbations is to 
liberate the inflaton from the task of generating 
the observed amplitude of perturbations. 

In (see also PHI) it was shown that, if cos¬ 
mological perturbations are produced during the 
late decay of a light scalar field whose value dur¬ 
ing inflation is of the order of mp, then the pri¬ 
mordial curvature perturbation already present 
from inflation, i.e., (),, is of the same order of 
magnitude as the perturbation produced by the 
late decay. This is true only if inflation is of slow- 
roll type, and fluctuations in the inflaton field are 
~ 77*, i.e., of the same amplitude as the fluctua¬ 
tions in the light field </. It is thus interesting to 
compare these two contributions in the context of 
the mass domination. 

The curvature perturbation /* from inflation 
is, at first order in the slow-roll parameters, /* = 
(l/V^)(77*/mp) [T], where e <C 1 is the first 
slow-roll parameter. We can replace this expres¬ 
sion in the definition of R, Eq. 


r -. 

3 /)* 


(39) 


Unless 0* is very small, 0* ^ fy/em-p, the pertur¬ 
bation generated by the mass domination repre¬ 


sents only a negligible correction to the perturba¬ 
tion originally present from inflation. However, 
we know that (/* cannot be too small: Indeed, 
Eq. 0 implies that (/*/mp > which com¬ 

bined with Eq. (PI yields 

R < (40) 

The stronger condition cni implies an even lower 
value of R. We conclude that if inflation is of 
slow-roll type, the mass domination mechanism 
does not provide a sufficient amplitude of fluctu¬ 
ations to liberate inflation. 

This conclusion however changes if the infla¬ 
ton does not go through a standard slow-rolling 
phase, as proposed in |3]. If the inflaton is not 
light during inflation, its fluctuation is suppressed 
m and the primordial curvature perturbation 
originated during inflation is negligible, /* ~ 0 . 
Only in this case, the inhomogeneous mass domi¬ 
nation mechanism represents a viable scenario for 
the generation of primordial cosmological pertur¬ 
bations. 

8. Conclusion 

We have considered the possibility that cosmo¬ 
logical perturbations are generated by the mass 
domination mechanism, as proposed in in, dur¬ 
ing a phase of domination of non-relativistic par¬ 
ticles whose mass is fluctuating in time and space, 
modulated by a light scalar field. By requir¬ 
ing that the scalar field remains light until the 
massive particles decay and does not dominate 
the universe with a second stage of inflation, we 
have shown that its value must be of the order 
of the Planck mass. A further condition must 
be imposed on the annihilation cross section of 
the massive particles. This must be weak enough 
- weaker than for standard model interactions - 
as to let the massive particles freeze out before 
being completely diluted when they become non- 
relativistic. 

We show that if these conditions are met, a cur¬ 
vature perturbation is produced. This is propor¬ 
tional to the mass fluctuation and the abundance 
of massive particles at their decay, as given by 
Eq. 115511 . and confirms E], in the limit where the 
massive particles completely dominate the uni- 



verse. Non-Gaussianities are inversely propor¬ 
tional to the abundance of massive particles at 
their decay, so that the latter can not be too 
small. 

Perturbations produced during inflation add to 
those generated by the mass domination. If in¬ 
flation if of slow-roll type, the former are much 
larger than the latter. Only for inflationary mod¬ 
els that violate the slow-roll conditions is the mass 
domination a successful mechanism for generating 
the observed cosmological perturbations. 
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